function variety_draws = Frechet_draws(T,theta,draws)

CDFvals = rand(draws,1);

variety_draws = T.*((-log(CDFvals)).^(-1/theta));
end
